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For a uniform strut whose moment of inertia is only 1/4 that at the mid- 
section of the above strut 
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Hence by tapering the strut to a moment of inertia at the tip equal to 1/4 
that at the central section, the strength has been reduced by 16 per cent. 
The variation in strength does not appear to be greatly affected by moderate 
tapering. The last formula, therefore, should cover most practical cases. 
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Two-dimensional Solutions of Poissons and Laplace's 

Equations, 

Part I. — By Leonard Bairstow, F.E.S., C.B.E. Part II. — By Arthur 
Berry, M.A., Fellow of King's College, Cambridge. 

(Received January 13, 1919.) 

In the course of a preliminary investigation of fluid motion when viscosity 
is taken into account, it was found that a solution of Laplace's equation was 
obtained in very general form. The solution gives exact analytical expres- 
sions where the conditions of the problem admit of exact solution, but the 
chief claim to notice arises from the extremely general character of the 
method if graphical and mechanical integration are employed. 

The problem has been approached from the engineering standpoint, and 
the absence of any proof of the convergency of the series developed will 
immediately be evident to the mathematician. The justification by trial in 
particular cases is all that is attempted at the present moment, but in those 
cases the result is wholly satisfactory. For the analytical solution of the 
problems I am indebted to Mr. Arthur Berry who has added a second part 
to this paper, and has generally acted as helpful critic. 

For single closed boundaries of any shape, or for a combination of separate 
boundaries, the problem of finding a function which satisfies either Poisson's 
or Laplace's equation, admits of direct attack for the region inside a single 
boundary or between boundaries if there be more than one. For the region 
outside a single boundary there is a further condition to be satisfied infinitely 
far from the boundary. 

The problems arising in calculating the torsion of solid or hollow cylinders, 
the flexure of beams, and the motion of an inviscid fluid round a cylinder 
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are covered by the new analysis. The functions developed appear to belong 
to the class called ''Greenes functions," and are special to the boundary or 
boundaries involved in any particular problem. 

The nearest approach to the degree of generality now attained is to be 
found in the application of soap films to problems of torsion and. flexure. 
An account of that method is given by Griffith and Taylor in a paper read 
before the Institution of Mechanical Engineers in December, 1917. A 
visible representation is given which may often be translated directly into 
an expression of the strain in a cylinder or beam. All problems which 
may be solved by the use of soap films may also be solved by the new 
analytical method, which is further capable of attacking problems where the 
soap film is unsuitable. 

The basic theorem used is given in Lamb's ' Hydrodynamics,' and in Love's 
' Elasticity ' as due to Green. In the language of hydrodynamics the theorem 
states that " any irrotational acyclic fluid motion can be reproduced by an 
appropriate choice of simple sources round the fixed boundaries." The more 
general statement given in the works referred to is not required in the present 
application. 

The strength of the sources corresponding with given boundaries is deter- 
mined by the solution of an integral equation. In the case of double 
boundaries simultaneous integral equations are solved. 

The general form of Poisson's equations in two dimensions is 

V^fi=/(^.?/)> (1) 

where /is a known function, and where yjri or d'>lrildn has a known value at 
every point on the fixed boundaries. It is known from the theory of attrac- 
tions that a particular integral of (1) is given by 

^2 = 2^ I J /(^> y) log T dxdy, (2) 

where r is the distance from the point at which ^2 is to be determined to 
the point {x, y). It is then permissible to con&ider the equation 

as the general case with the boundary values of either i|r or dyfr/dn given. 
Equation (3) can be solved in terms of sources on the boundaries if the 
motion is acyclic, and if an extension to cylic motions is required, this can be 
attained by the introduction of artificial barriers as has been pointed out by 
Lamb. The extension will not receive further notice in this paper. 

To find the source distribution round a single boundary, fig. 1, it is 
necessary to build up the efifect of a number of elementary sources. At Q, 
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the elementary source d% adds to the stream function -^p an amount 
1 



Jtt 



dd%, and hence 



-^p = 2^ I &^l^. 



(4) 



where the integral is taken all round the boundary. % is a function of 



FfS.L 




(fig. 1), and is indeterminate to the extent of a constant. Change the inde- 
pendent variable to 6 by partial integration, so that 



"^p = 



27rL 



e% 



= 2Tr 



= 
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Xdd. 



(5) 



If X be not discontinuous at the limits, the expression between limits leads 
to a constant, and without loss of generality we may write 



^1 



27rJ 



Xde. 



(6) 



Now consider equation (6) when the point P lies on a fixed boundary. The 
integral can be divided into two parts, one of which is equal to — tt^p and 



p0] + 7r 



the other is Xdd, where @i is the angle of the tangent at P. The positive 



J0 



direction is to the right when standing on the boundary facing into the 
enclosed region. The particular feature of the integral here referred to is 
seen by considering the limiting conditions as P aijproaches the boundary 
from outside. Equation (6) becomes 



1 r®+7r 
yfr — 1 % — -— I Xdd, 



(7) 







Had the point P been inside the boundary the expression corresponding 
with (7) would have been 
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2'^P 



ZTT J 



(0 + 7r 
xde. 



(8) 
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Equations (7) and (8) are integral equations, of which the solutions are 



1 

IT 



"© + 7r J[ /'0+7rr0 + 7r 



and 



ix, = ,|r,+- irde+—A \fdedex+..., (9) 

l%p=: ~-fp+i ylrdO-W fd0dei+,.,, (10) 

ttJ® tt j@ J© 

the two series being extended indefinitely. Each integral is taken for Q 
travelling round the boundary, but the radius vector traversing only an 
angle tt. The double integral signifies the operation repeated on the result 
of the first. 

If the series is convergent the solution is easily checked by substitution 
of (9) and (10) in the original equations (7) and (8). 

An examination of the terms of (9) will show that the first integral is an 
average value of sjr on the boundary. The second integral is an average of 
the same kind on the results of the first integration, the base on which the 
average is taken being the angular position of all points Q relative to P 
(fig. 1). The result in the particular cases tried is that the successive 
integrals tend to become independent of the position of P but are not equal 
to zero. 

If, as an example, we assume that -il^p is constant and equal to C, 
equation (9) becomes 

l%p=: C + C4-C+.,., (11) 

and the series is not convergent. There is, however, no objection to the 
expansion of 'v//^p—C in the solution of equation (3), and the difficulty shown 
in (11) is then removed. 'No physical explanation of the meaning of C in 
the hydrodynamical problem is given, but the existence of a constant 
residue C does not introduce any difficulties in numerical computation, as it 
is automatically evaluated in calculating the successive integrals on the 
assumption that it exists. 

Fii'st Uxmiiple of the Use of (8). 

If an inviscid fluid be still at infinity, and a circular cylinder be moving 
through it with velocity V, we have the boundary condition 

fp = YaBm@, (12) 

where a is the radius of the cylinder and @ is measured from the centre of 
the cylinder, and is zero in the direction of motion. 
Applying equation (9) to (12) leads to 

|Xp === Vasine + ^fsin@6^^+™^ffsin(H)^^^0x+...^ ' (13) 



Solutions of Poissoris and Laplace s Equations, 46 1 

and since dd = ^d&y the integrals can be evaluated and found to be zero 

separately. We then find 

Xp= 2Vasin@. (14) 

The value of sfr at any point not on the boundary is now found from the 
equation 

^= ^^Lin®de. (15) 

This integration presents no particular difficulties and leads to 

y}r = l±.^inei, (16) 

r 

where r and Oi are polar co-ordinates measured from the centre of the 
cylinder and the direction of motion. This solution is of course well known. 

In the general case of a boundary of any shape, the integrals in (9) 
and (10) cannot be obtained analytically. If graphical and mechanical 
integration are introduced they are always calculable, the degree of accuracy 
being fixed by the care bestowed on the arrangement of the work. As the 
examples now to be given are illustrative only, no attempt has been made to 
achieve more accuracy than is sufficient as an effective check on the con- 
vergence of the series. 

The remaining examples refer to the torsion of cylinders, and the equation 

to be solved is 

V't = 0, (17) 

with the boundary relation 

'^6 — 2 (.^^ + y^) = const. (18) 

Second Example : The Torsion of an Elliptic Cylinder, 

The origin of co-ordinates was taken at the centre of the ellipse with the 
axis of X along the major axis. The lengths of the semi-axes were 1 
and 0*5. Points on a quadrant were marked 0, 1, 2, 3, ... 10 to denote 
position; point was at the end of the minor axis and 10 at the end of the 
major axis, the intermediate points having abscissae equal to one-tenth of 
their marking. 

The numerical values of the several terms of equation (10) are shown in 
Table I, and it will be noticed that the second integral has attained a value 
which is little different from 0*317 at any point of the boundary. All 
further integrals of the series would have this same constant value, and the 
series will be convergent if the constant of equation (18) is made equal 
to — 0*317. Using this value leads to Table II, which shows how very 
rapidly the series converges in this instance. The value of %, i.e., the 
integral of the strengths of the sources on the elliptical boundary, is given in 
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the fifth column of Table 11, and, for comparison, a column is given of the 
values computed analytically by Mr. Arthur Berry (see Part II). The 
agreement is complete within the accuracy of the graphical calculation. 

Third Example : The Torsion of a Cylinder of which the Section is that of an 

Airscrew Blade Section. 
The section of the cylinder is shown in fig. 2. The origin has been taken 
near to the centre of area, and the points 1-11 were chosen arbitrarily. 
The figure gives only the shape, the values of |r^ at each of the points being 
given in Table III, column 2. 




After three integrations, it appeared that the limiting constant value of 
the integrals was approximately 45* 1^ and, using this value for the boundary 
constant of (18), the calculated value of successive terms of equation (10) are 



given in Table III. 



Table I.-— -Elliptic Cylinder. 



Location of 
point. 


i.e., i|/(j — const. 


-- ) (i|/^; — const.) d9. 
TV «' 

0-296 


—^ If (+(!' — ^oiist.) dd dBi. 





-125 


0-316 


1 


-129 


0-296 


0-316 


^ 2 


0-140 


0-298 


0-316 


3 


0-159 


0-299 


0-317 


4 


-185 


-301 


0-316 


5 


-219 


0-304 


0-316 


6 


0-260 


0-307 


0-316 


7 


0*309 


0-313 


0-318 


8 


0-365 


-321 


0-318 


9 


0-429 


1 -330 


0-318 


10 


0-500 


i -336 

i 


-318 
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Table 11.™ 


"Elliptic G^ 


binder. 




Location 

of 

point. 


(1) 


(2) 




^2[-(l) + (2)-(3)]. 


X as deduced 

analytically 

(see Part II). 





-0 'im 


-0-021 


-0-001 


+ 0-344 


+ -345 


1 


-0-188 


-0 -021 


-0-001 


+ 0-336 


+ -338 


2 


-0-177 


-0-019 


-0 -001 


+ 0-318 


+ 0-318 


3 


- 0-158 


-0-018 


-0-000 


+ -280 


+ -284 


4 


-0-132 


-0-016 


-0-001 


+ 0-234 


+ 0-2'37 


5 


-0-098 


-0-013 


-0-001 


+ -172 


+ 0-176 


6 


-0-057 


-0-010 


-0 -001 


+ 0-096 


+ 0-102 


7 


-0-010 


-0-004 


+ 0-001 


+ -010 


+ 0-014 


8 


+ -048 


+ 0-004 


+ -001 


- -090 


-0-087 


9 


+ 0-112 


+ 0-013 


+ -001 


-0-200 


-0 -202 


10 


+ 0-183 


+ 0-019 


+ -001 


-0-330 


-0-330* 



^ The constant of integration in the analytical solution was chosen to make this figure agree 
with that calculated graphically. 

Table III.— Cylinder of which the Section is that of an Airscrew Blade 

Section. 



Location 
of point. 


J(a^H5^2). 


^b i(ajH/)-45'l. 


i {y^bdB. 


i f Uz, dd ddi. 


-^ Ul\pi,ded0id02. 


1 


89-7 


44-6 


24-6 


13-1 


7-3 


2 


52-2 


+ 7-1 


+ 5-8 


4-0 


2-4 


3 


9-6 


-35-6 


-10-4 


-3-0 


-0-8 


4 


23-1 


-22-0 


- 8-5 


-2-9 


-1-7 


5 


45-6 


+ 0-5 


- 5-1 


-1-7 


-1-9 


6 


86-7 


41-6 


+ 10-4 


+ 3-5 


-0-2 


7 


79-9 


+ 34-8 


10 -3 


+ 1-6 


0-0 


8 


28*3 


-16*8 


- 4-5 


-2-9 


-1-2 


9 


9-7 


-35-4 


- 9-9 


-3-5 


-1-3 


10 


24-1 


-M-O 


- 6-5 


-0-9 


0-0 


11 


68-0 


+ 23-9 


+ 12-7 


+ 7-7 


+ 4-7 



The convergence is not nearly so rapid as for the ellipse. The trend of the 
terms is clear, and it would .be easy to guess an extension, and so approxi- 
mate to the value of X very rapidly. This is always permissible, since the 
repeated processes do not presume correctness of the previous work. The 
value of y^if, calculated from the approximate value of X, will have small 
errors in it, which will lead to a new problem with new boundary conditions. 
The importance of this possibility of approximate calculations is very great, 
since the extra care necessary for high accuracy can be introduced at a late 
stage of the computation. The calculations were made by Miss A. D. 
Betts, B.Sc. 

To the accuracy given in Tables I, II, and III, the graphical process is 
particularly simple, and each of the examples 2 and 3 can be solved in a 
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few hours by a mechanical draughtsman. Once X has been determined with 
the desired degree of accuracy, it is possible to calculate i|r at any point of 
the section. Usually the stresses are required, i,e,, d^/dx and d'^jdy, and 
these may be estimated independently of y^. The mechanical operations for 
finding dXjdoo and o^jdy are not so simple as for -^^ but they present no 
particular difficulty. In the hydrodynamical simile used throughout this 
paper, the calculation of stress is equivalent to the calculation of the 
velocity of the inviscid fluid. 

Solution of Mqtmtion (3) taken d^jdn on the Boundary is Knoivn. 

If dy^Jdn represent variation normal to the boundary, the problem can be 
transformed by the use of Green's theorem so as to reduce to that already 
solved. If -^ be any harmonic function, Green's theorem leads to the 
equation 

27ri|rp== -~|2xiog^^5 4.L^^^^^^ (19) 

where 6i is the angle between the line PQ (fig. 1) and some fixed datum, 
and cU is an element of length of the boundary. With the point P on the 
fixed boundary, equation (19) becomes 

^ +-1 ^ItiycWi = — I ^logrds, (20) 

^ TTJe TTicn ^ 

the ambiguity of sign corresponding with solutions for the regions external 
to and internal to the fixed boundary If we now refer to equation (9), it 
will be seen that it may be written as 

ip^p = /^p + 1 U^fZ^^il+i \ ^^^+-{^tde^de2d0'i..., (21) 

Since the right-hand side of (20) is known by hypothesis, a combination 
of (20) and (21) leads to a determination of %, and therefore to a solution of 
equation (3), with dyjr/dn known on the boundary. 

Extension of the Analysis to Double Boundaries, 

Double boundaries can be treated by an extension of equations (9) and (10). 
The computation is of the same nature as before, but the number of integrals 
to be evaluated is greater. In developing the formulae, %o and %i will be 
used to denote the value of X on the outer and inner boundaries respectively. 
It is necessary to have a notation which will indicate that the points P and 

Q are on one or other of the boundaries. has been used to mean that Q is 
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oil the outer boundary, and that Q is on the inner boundary. denotes 

that P is on the outer boundary, whilst refers to P on the inner boundary. 

The notation does not lead to confusion with the limits of integration, 
since Q is always taken to traverse the boundary on which it lies. 
With the notation explained above, 



^fr, = ^[\ode+J^\Xide, 



27r 



2 



TT 



(22) 



If P be on the outer boundary, (22) becomes 



TT 1^ 27r 



(23) 



whilst for P on the inner boundary (22) becomes 

-.y^i = l-{\cie+~\Xide. 



(24) 



A usual boundary condition will be that ^^o and y^i are known functions 
except for additive constants. Equations (23) and (24) then constitute 
a pair of simultaneous integral equations from which to determine Xo and X^-. 
They are solved as follows : — 



Write X10 + X20 for X^ and substitute in (23) and (24) to get 

~^, = A \\od6 + ^ [\od9 + ^[x,de, 



and 



1 fo '] ro 



1 f\ 

X2od6 + -^ — ' 
i 27rJ. 



Xidu. 



(25) 



(26) 



Now define Xio by the relation 



■^0 = , 



2ir 



XiodO, 



(27) 



i.e., by an equation which is identical in form with equation (6), and which 
can be solved for %io, as previously indicated. Two equations are then left 
instead of (25) and (26), and are 



i^i 



'>, 



^ X2odd + 

J 

1 ro 1 ro 

X,„dd =: 



Xidd, 



IT 



■2tt 



X2odd + :^\Xid0. 



(28) 



(29) 



Since Xi^ has been determined by (27), we may look on the left-hand side 
of (29) as known, and write equation (29) as 



fu = T^ 



1 



TT 



X20 djQ -^ 



2 



TT 



XidO, 



(30) 
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an equation of exactly the same form as (24), and therefore suitable for 
similar treatment to (23), ie., write %it-f-%2t for %», and substitute in 
equations (28) and (30) to get 



= 



ro 



rt 



X^odd+\Xiode+ X^icte (31) 

Jo Jo Jo 



d -fii = 7r{^^2od0~{~~~{\dd+,^^^ (32) 

^'^Ji ^TTJ,- ^TTJ^. 

Define Xn by the relation 



1^ 



-^it — :y;^ I XiidO, (33) 



This completes a cycle of operations which can be repeated indefinitely. 
To see the repetition most symmetrically, write 

1 r 

"-i|r2o for Tp- XiicW, (34) 

when (31) becomes 

- f2o - -^hode + ~;^\x2id0, (35) 

whilst equation (32) becomes 

=:= ^hode + ^ fc^^. (36) 

These equations are essentially of the same form as (23) and (24). The 
zero term of (36) implies a considerable reduction in the labour of computa- 
tion and follows from the combination of the two terms of equation (29) into 
a single one in (30). Had this not been done equation (36) would have had 
still more exactly the form of (24). 

If X20 = Xso + X^o, and %2i = ^3i-f^4i, (37) 

and so on in the repetitions of the process, then 

Xq = Xio -f X20 -f- X^o + . . . "^ 

r- (38) 

The physical interpretation of the analysis in terms of fluid motion is not 
difficult and is given below :— 

(1) Find the distribution of simple sources on the outer boundary which 
would give -v/r^ its correct value if the inner boundary did not exist 
(equation 27). 

(2) Find the value of i|r on the inner boundary due to the sources found 
by (1). Add this value of yfr to the value of yj/i on the inner boundary (left- 
hand side of 29). 
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(3) Find the source system on the inner boundary which gives to 'y^ on 
the inner boundary the combined value of '^i given by (2), ignoring the 
existence of the outer boundary (equation 33). 

(4) Find the value of ^ on the outer boundary due to the sources of (3). 

(5) Find the source system on the outer boundary which eliminates the 
value of i/r due to the sources on the inner boundary, ignoring the existence 
of the inner boundary in the calculation. This is really a repetition of (1). 
Continue the processes until the various integrals become sensibly constant 
at each stage, and eliminate the constants as shown after equation (11). 

(6) Add the several solutions on the outer and inner boundaries respectively 
to get the final values for the source system on the two boundaries. 

Fourth Example, 

The torsion of a hollow circular cylinder is considered, the inner cylinder 
being circular but eccentric to the exterior surface. The radius of the outer 
cylinder was 5 and of the inner 2*5, the distance between their centres being 
unity. 

The origin of co-ordinates has been taken as the centre of the outer circle 
and the location of the points referred to in Table IV is given in terms of 
the angles measured from the centres of the respective circles in the section 
of the cylinder. . The points and 9 are those of greatest separation on the 
two circles and points 8 and 17 those of nearest approach on the line through 
the two centres. The angular separation of points on either boundary is tt/S. 
The symmetry of the section about the line joining the centres of the circles 
makes it unnecessary to calculate integrals on both halves. 

The equation to be solved is 

with -xIto— |(a?2+^2) = const. > (39) 

■^i— 'K^^H-?/^) = const. ^ 

Griffith and Taylor in solving (39) by the soap film method found it 
necessary to introduce the condition that the conjugate function ^ is to be 
single valued. The hydrodynamical analogy is that the motion between the 
boundaries shall be acyclic and this condition is automatically fulfilled in the 
present analysis. 

Since the centre of the outer circle has been taken as origin, '^o is constant 
and Xio = 0. On the inner boundary 

'fi rr I {r^ ^cPj^2 Ti d cos )8) -I- const., ' (40) 

d being the distance between the centres of the circles of the section and /3 
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the angle measured from the line joining the centres, using the centre of the 
inner circle to fix the axis of rotation of the radius vector. 

Table IV.— Hollow Shaft. Section made of two Eccentric Circles. 



Location 

of poinfc, 

outer circle. 




1 
2 
3 
4 
6 

7 
8 



X 



10' 















Location 

of point, 

inner circle. 



9 
10 
11 
12 
13 
14 
15 
16 
17 



^1,-. 



^3o« 



X, 



Vif 



>^5.- 



-2 -04 


-2-00 


-1-72 


-1-24 


-0-48 


+ 0-48 


+ 1-58 


+ 2-64 


+ 2-98 



(-0-50) 

(-0-48) 

(-0-40) 

(-0-32) 

(-0-12) 

( + 0-10) 

(0 -42) 

(0 -68) 

(0 -74) 



^51. 



-5-00 


-1-18 


-4-66 


-1-12 


-3-60 


-0-92 


-2-00 


-0-58 





-0-12 


+ 2-00 


+ 0*42 


3-60 


+ 0-92 


4-66 


1-30 


5-00 


1-44 



(-0-28) 

(-0-28) 

(-0-24) 

(-~0-14) 

(-0-04) 

(0 -lO) 

(0 -24) 

(0 -32) 

(0 -36) 



X. 



<?• 



-9o« 



(--0-12) 
(-0-12) 
(-0-10) 
(-0-08) 



( 



0-04) 
(0 -04) 
(0 -10) 
(0 -16) 
(0 -18) 



(-0-02) 

(-0-02) 

(-0-02) 

(-0-02) 

(0 -00) 

(0 -00) 

(0 -01) 

(0 -04) 

(0 -04) 



^<7. 



from 

Xg. 



X' 



-;:z:2::;-;::;.-_:-^::i^ 


1 


^7/. 




(-0-06) 


i 


(-0-06) 


j 


(-0-06) 




(-0-04) 




(0 -00) 


i 


( + 0-02) 




(0 -06) 




(0 -08) 




(0 -08) 


1 
1 



-2-68 


0-08 


-0-16 


-2-62 


0-08 


-0-16 


-2-24 


O-08 


-0-16 


-1-66 


0-06 


-0-12 


~0-64 


0-02 


-0-04 


0-62 


-0-03 


+ 0-06 


2-12 


-0-09 


0-18 


3-52 


--0-15 


0-30 


3-94 


"0-18 


0-36 



Xq as deduced 
10. analytically 
(see Part II). 



-2 -82 



X,. 



-6 
-6 

-4 

-2 

-0 

+ 2 

4 

6 

6 





t/^z + k^ 




-12-5 




from 




X,. 


52 


-8 -60 


12 


-8-58 


82 


— 8-57 


72 


-8 -57 


16 


-8 -61 


54 


-8 -70 


82 


-8-61 


36 


~8-56 


88 


-8-50 



-2 
-2 
-1 

-0 

+ 
2 
3 

4 



•70 
'34 
•64 
•64* 

-m 

•34 

•82 
•46 



Xi as deduced 

analytically 

(see Part II). 




* Made to agree with Xo as calculated graphically by choice of constant of integration, 
f Made to agree with X,; as calculated graphically by choice of constant of integration. 



Analytically, it is found that 



l%ii = nd cos ^, 



(41) 



%3e and Xgi were determined by graphical and mechanical integration and 
are given in Table IV. It was clear at this stage that the series for Xo and 
%t were both convergent with fair rapidity and the main object of the 
illustration was achieved. In order to get a rough approximation to %o and 
Xi remaining terms were estimated by observing that X^i is roughly one 
quarter of Xn and that as a consequence X^^ will be one quarter of X^i and so 
on. The values of Xo andXt so obtained are given in Table IV. An estimate 
of the residual error was then made by calculating yjto and -v/r^-f -|^^ — 12*5 
both of which should be constant. In both cases departures from constancy 
were observed and X\o was calculated as the next approximation in a new 
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series by doubling '>^o and changing sign.* This gives a good indication of the 
errors in the first value of Xq- 

Mr. Arthur Berry has calculated %o and Xi analytically, and the values are 
given in the last column of Table IV. The approach to these values by 
%o and %i is not so good as for the elliptic cylinder, but the estimate of X'u 
shows that this approach could be improved rapidly by further calculation if 
any use were to be made of the results. 

Conclusion, 

From an engineer's standpoint it is thought that the theorem is reasonably 
established and may be used with confidence in new problems, since any 
answer obtained constitutes the necessary proof of convergency. It is hoped 
that mathematicians will complete the work on their side. It appears from a 
further theorem given in Lamb's ' Hydrodynamics ' that the equation 

may possibly be capable of solution along the same lines as the present 
solution of Poisson's equation, the elementary source being changed for the 
new equation. 

Part II — By Arthur Berry, M.A., Fellow of King's College, Cambridge. 

§ 1. Introduction, 

At Mr. Bairstow's request I have carried out the analytical work corre- 
sponding to his method of solution for two cases of the torsion problem in 
which there is a known analytical solution, viz., the elliptic cylinder and the 
hollow beam bounded by non-coaxial circular cylinders. 

We interpret the torsion function as the stream function of two-dimensional 
irrotational fluid motion, defined by <!>-{- i-yjr, within a region corresponding to 
the section of the beam. 

In the case of a simply connected region we have to determine for the 
external space a continuous acyclic fluid motion defined by a velocity 
potential ^' such that ^'vanishes at infinity and agrees with ^, save as to an 
additive constant, on the boundary. 

The integral ^^l^lP^ P^) ^'> (1) 

27rJ \dv ov I ^ ^ 

taken round the boundary of the region {0 being the angle which the radius 

^ For a circular boundary the source system is equal to twice the boundary value of 
^ except for an additive constant, no matter what the source system may be. 
VOL. XCV.— A. 2 P 
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vector from a variable point, xy, on the boundary to a given point, x'y\ in the 
interior of the region, makes with an arbitrary initial line) is the stream 
function at x\/ due to a distribution of sources on the boundary of strength 
d<p/dv-\-d(f>^ /dv' per unit length, and is known to agree with '>|r save for an 
additive constant. The sense of integration is the same as that in which & is 
measured. 

In Mr. Bairstow's notation -— = ;=^-|-^ . 

as OP ov 

For a doubly connected region we have similarly to determine (p'o for the 
external space, so that </)'o vanishes at infinity, and agrees with (j> save for an 
additive constant on the outer boundary, and also ^\- for the internal space 
so that (^'i agrees with cj> save for an additive constant on the inner boundary,, 
the fluid motions defined by these functions being continuous and acyclic. 
The stream function at any point of the interior of the doubly connected 
region is then given by 






TT 



^&-^--k f"©-'^^*. «■ 



where the two integrals are taken round the outer and inner boundaries 
respectively in the sense in which increases. In Mr. Bairstow's notation 

ds av dv' ' ds dv 3/ 

In both cases, since by the conditions of the problem the sum of the 
strengths of the sources on either boundary is zero, the integrals are inde- 
pendent of the initial line from which 6 is measured, but as ^ is a cyclic 
function they are not independent of the point on the boundary at which 
the integration begins, so that each involves effectively an arbitrary constant 
of integration depending on the choice of the lower limit of the integral. 
The stream function calculated in this way will, in general, differ by a 
constant from the function -v/r determined by other methods, e.g., determined 
to satisfy the condition 'f = ^{x^ + f) on one of the boundaries. This is a 
partial explanation of the fact that in Mr. Bairstow's work it is necessary to 
add a constant in order to ensure the convergence of the series which 
he uses. 

§ 2. The Elliptic Gylinder. 

The semi-axes being a, 6, we can take 

.|r = X(a^2_^2)^ ^=^2\xy, (3) 

where X = i(a^-&^)/(aHn 
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since ^, i|r, are conjugate harmonic functions satisfying the continuity con- 

ditions, and 

i|r = 1 (a?^ + if) --a%^l{a^ + W) 
on the boundary. 

Introducing elliptic co-ordinates defined by 

^ = c cosh f cos 7;, y =z G sinh f sin ^;, 

where c^ = a^—5^ and f = ^0 on the boundary, we have 

^ = — I c2^ sinh 2 f sin 2^, 

i|r = I c^X (cosh 2 f cos 2'^ + 1). 

We can now satisfy the conditions of the problem by choosing 

^;.= -ic2xg-2^'-2^sin29;sinh2^o, (4) 

since this is a one-valued and continuous harmonic function which vanishes at 
infinity and agrees with <l> on the boundary. 
We now have 



as \dv ov 










so that 



^X (cosh 2 ^0 + sinh 2 f 0) sin 2 9; % = | X (a + &)^ sin 2 <?; c?-?;, 
X = — ,, — -.:(c^-f6)2(l— cos2'^), 



the constant of integration being chosen so that X vanishes with ^. 
In the numerical case chosen by Mr. Bairstow, in which a = 1, 5 = |-, 

27 



(5) 



% 



80 



(1— cos 217), 



and the numerical values at the eleven points chosen by him are given by 
the Table : — * 



Point. 



cos 97. 


-6750 





1 


0-1 


-6682 


2 . 


0-2 


0-6480 


3 


0-3 


"6142 


4 


0-4 


-5670 





0'5 


-5062 


6 


0-6 


-4320 


7 


0-7 


0-3442 


8 


0-8 


-2430 


9 


0-9 


0*1282 


10 


1-0 






* I am indebted to Miss E. B. Lang, M.A., for the calculation o£ 

this Table and of those in §3. 

2 p 2 
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In this simple case it is easy to carry out analytically the successive steps 
of Mr. Bairstow's approximation. 

Take i/r = -| {p(? + y^) on the boundary. 

If Xoyo is any point on the boundary which is fixed during the integration, 
we can put for any other point on the boundary 

X = Xo-\-r cos 0, y z=: y^-^-T siu 6, 

r, 6 being polar co-ordinates referred to Xo, yo as origin. 

cos^ 6 . sin^ 6 



mi o {^^0 cos 6 , Vo sin 0~1 / 



a' 



¥ 



and 



1 

IT 



fcW 



1 TTT + a 

{Xo^ + yo^-{-2 XqT cos 6-\-2yoT sin O-^i^} d6, 



TV 



a 



where a is the initial value of 6 when x, y coincides with Xo, yo- 



ITow 



T^cW ~ twice area of the ellipse = 27rah 



r cos Odd 



Cir 



[xoQ.0^^6 sin cos 01 //cos^0 , sin^0 



ZiXi 



a' 



cos^ 9 sin^ 6\ ^^ 



- 



+ 



\ a^ 



W 



de 



•TV 



cos^0 



W 



— 27rXob/(a-\-h), 



and, similarly, 

1 '' 

Hence — 

TT 



r sin 6^0 = — 2 7ryoCc/(a + &). 



-^^0 = i(xo^^yo^)-2xo%/{a + h)''-2yohc/(a + b) + ah 

since x^ja^-k-y^jW = 1. 

If now we drop the suffixes, we can write this integral for brevity in the 

form 

Ii = fil{x^-]ry^)-\-v, 

where ^ = {a — hf/{a + by, 'v =■ ah {a? + 1^)1 {a + Vf. 

Eepeating the process, we have 



L = 



IT 



IT 



I3 = \\l2dd = /.2Ii + (/i+l)^ = |-/x3(«'+?y')+(M'+M+l)f' 



Since < /i < 1, I„ converges to a constant limit 

I„ =z./(l-/.)=i(aHn 
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Since this limit is not zero, Mr. Bairstow*s series — i|r + Ii — Ig-f... does not 
converge. As, however, any constant addition to i|r produces the same 
addition to each of Ii, I2, ..., we have only to take for the torsion function 
^ = -xlr— loo = '\{r—l{a^-{-h^), and we then obtain 

I2 = ^Ii = /J?^j?, 



In = fM'''\{r, 

and the series then converges to the sum 

which agrees save as to an additive constant with the form (5) previously 
obtained for %. 

We thus see that the subtractive constant which appears in Mr. Bair- 
stow's arithmetical work is i(c(P + b% or in his numerical case f^g- = 0*3125, 
agreeing closely with his result. 

§ 3. The HollpiD Beam lounded hy tivo Circular Cylinders. 

The solution of the torsion problem has been given by Macdonald* in the 
form 

I n 9 ^ / N« ^""^^ coth B sinh ^^ (w — a) 4- ^~^* coth a sinh n (B—v) 

^ = 2c^z( — r . T .^ r -^^ — -GosnP, 

I smh^(p — a) ^ 

where x + iy = c tan |(f + irj), 

and rj = a represents the outer boundary, rj = ^ the inner. 

Writing this for brevity, 

00 
'\/r = S {A„sinh?^(7; — a) + B^sinh^(/3— ^7)} cos^i^, 



1 

CO 



we have ^ = S {A^cosh?i('7?— a) — B^cosh^(/3— ^)} sin^^f. 

We can now take for the external region, ?; < a, 

00 

^/ = S6'^(^-'*){A^~B^cosh7^(/3-a)}sin^|^, 
1 

since this is a continuous acyclic harmonic function, which vanishes at 

infinity (^ + ii] = if) and agrees with ^ when rj = a. 

For the internal region, -?; > /3, we can take, similarly, 

<f)i' = i e-'^(^~^){A^cosh^(/3-a)~B4 sin^i^, 
1 

* ' Cambridge Phil. Soc. Proc./ vol. 8 (1893) ; Love's "Elasticity," § 222, Ex. (iii). 
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since this is a continuous acyclic harmonic function throughout the internal 
region and agrees with ^ when '^ = yS. 
We now have on the outer boundary 

as \ov OP I \07} dr) /,,=« 

= %n{Bn sinh n O — a) — A^ 4- B^ cosh n (/3 — ce)} sin n ^d^ 
1 

CO 

= 2 ^ { - A,, + B,,e^<^-«>} sin n^d^, 
1 

whence, choosing the constant of integration so that Xo vanishes at the point 
1^ = 0, where the two circles are closest together, 



CO 



^0 = 2 {-A„ + B„e»(^-«)}(l-coswf) 



1 

00 



o2x</ w »«,6'*^^"^^cotha-6-'^C^-«)coth^,, ^, ,^, 

1 ^ ^ smh.n{/3—ot) ^ ^^ ^^ 



Similarly^ on the inner boundary 



^*+^v^ 



3^7 d7] 
2^^ {—An sinh^iO— a)— A«cosh^i(/3--a) + B„} sin n^d^ 



1 

00 



The torsion function which we have used corresponds to a twist about the 

origin of co-ordinates chosen, which is the point of intersection of the radical 

axis and line of centres, and differs by a constant from l{x^-{-y^) on either 

boundary. Mr. Bairstow's torsion function corresponds to a twist about the 

centre of the outer circle (^ = 0, 2/ = ccotha), and differs by a constant from 

KX^ + Y^) on either boundary, where X, Y refer to the centre of the outer 

circle as origin. In order to obtain Mr. Bairstow's torsion function, we must 

therefore add to our ^^r the term — ccoth ay. 

Taking 

^|r = — c coth ay, (j) z= ^c coth ax, = — ccotli aX, 

we can take cj}o == — ccoth a :^^ = — c coth aTo^ cos 0/E, 

where To (= c cosech a) is the radius of the outer circle, and K© are polar 
co-ordinates corresponding to XY. 
We then have 

dX^ = ( -"*?r^ + -^ )d$ = c coth a cos ©ndM, 
\ oK oK / 

X^ = 2 c coth a To (sin @ + 1), (8) 
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where, as before, the constant of integration has been chosen so as to make 
%o vanish when © = — 7r/2 or f = 0. 

Since (^^ = <^, Xi = 0. 

For Mr. Bairstow's problem we have therefore to take for the inner circle 
%i, as given by (7), and for the outer circle the sum of the two functions %« 
given by (6) and (8), or denoting these for distinction by XJ and Xo\ we take 

%o = %/-!-%/'. (9) 

The points on . the outer circle numbered 0, 1, ... 8, correspond to 
€) = IT 1 2, 3/8, ..., — 7r/2, the corresponding values of f being given by 

cos f = —(1-1- cosh a sin 0)/(cosh a 4-,sin @), 
sin ^ == sinh a cos 0/(cosh a + sin @). 

For the inner circle we have corresponding formulae with ^ instead of a. 
In Mr. Bairstow's numerical case we have 

To == cQOBQohcc = 5, n = ccosech/3 = 2*5, . D = ccotha—ccoth/3 = 1, 

whence a == 1-3024 /3 = 1-9396, e = 8-5156. 

The values of X^, Xi are then given by the tables. 



Outer Circle. 



Point, 


©. 


■ I 


^0 . 




Ag, = A^ 4* ^0 • 





90° 


180° 


-204-78 


197 -50 


-7-28 


1 


67° 30' 


167° 1' 


-197-16 


190 -00 


-7-16 


2 


45° 


153° 20' 


-175-40 


168-60 


-6-80 


3 


22° 30° 


138° 8' 


-142-66 


136 -56 


-6-10 


4 


0'' 


120° 25' 


-103-86 


98-76 


-5-10 


5 


-22° 30° 


98° 50' 


-64-76 


60-96 


-3-80 


6 


-45° 


71° 46' 


-31-04 


28-92 


-2-12 


7 


-67° 80' 


38° 19' 


-8-16 


7-52 


-0-64 


8 


-90° 


0° 












Inner Circle. 



Point. 


1. 


Aj'. 


9 


180° 


-13-818 


10 


163° 3' 


- 13 -342 


11 


145° 83' 


-11-976 


12 


126° 51' 


-9 -884 


13 


106° 22' 


-7-332 


14 


83° 30' 


-4-662 


16 


57° .64' 


-2*278 


16 


29° 45' 


-0-606 


17 


0° 






